arXiv: 1502.07279v4 [quant-ph] 27 Jul 2015 


Fully QED/relativistic theory of light pressure on free electrons by isotropic radiation 


A. E. Kaplan 

Electr. & Computer Enqineerinq Dept., Johns Hopkins University, Baltimore , MD , 21218, USA*\ 

(Dated: July 28, 2015) 

A relativistic/QED theory of light pressure on electrons by an isotropic, in particular blackbody 
radiation predicts thermalization rates of free electrons over entire span of energies available in 
the lab and the nature. The calculations based on the QED Klein-Nishina theory of electron- 
photon scattering and relativistic Fokker-Planck equation, show that the transition from classical 
(Thompson) to QED (Compton) thermalization determined by the product of electron energy and 
radiation temperature, is reachable under conditions for controlled nuclear fusion, and predict large 
acceleration of electron thermalization in the Compton domain and strong damping of plasma 
oscillations at the temperatures near plasma nuclear fusion. 


I. INTRODUCTION 

Beginning with Max Planck discoveries [1], one of the 
fundamental issues in optics, electrodynamics, thermody¬ 
namics, atomic physics and quantum mechanics, is how 
a radiation, in particular blackbody radiation, imposes 
an equilibrium in the material system by either heating 
it up or cooling down. This process is facilitated by a so 
called light pressure [2-4] on charged particles [5], most 
of all the lightest ones - electrons, either (quasi)free as in 
plasma, of bound as in atoms or ions, in which case the 
electrons pass the light pressure on to atoms. The advent 
of lasers allowed the development of highly controlled and 
engineered non-thermal radiation environment, such as 
coherent laser light with its frequency tuned near atomic 
resonances, and use them for the cooling of atoms by 
red-shifted laser [6-9] (and coherent motional excitation 
of atoms by blue-shifted one [10,11]). A related pon- 
deromotive, or field-gradient force [12-14], manifested e. 
g. via Kapitza-Dirac effect [15-18], has also been used 
in laser trapping of atoms [6-9] and macro-particles [19], 
high-field ionization of atoms [19-21], etc. It can be very 
sensitive to relativistic effects, which under certain con¬ 
ditions may result in a chaotic motion [22,23] or even 
reverse the sign of that force in a strong field [24,25]. 

In the case of blackbody radiation acting upon free 
electrons, the situation is rather straightforward and fun¬ 
damental: non-resonant light pressure plays the role of 
the equilibrium ” enforcer” by either energizing slow elec¬ 
trons or damping the fast ones. The main issue here is 
how fast the equilibrium/thermalization can be reached. 
The relaxation time, t r i x , of that process could vary by 
many orders of magnitude depending of the temperature 
T and initial energy/momentum of electrons; in classi¬ 
cal domain t r i x oc T -4 . At the temperatures near ab¬ 
solute zero (e. g. in the so called relic radiations, or 
Cosmic Microwave Background, CMB [26-28]) the equi¬ 
librium is essentially unreachable (it is too long even at 
T ~ 10 3 A', see below), whereas in high-T environments, 
e. g. controlled nuclear fusion, nuclear explosions, and 
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star cores, it can be reached faster than in attoseconds. 
We found that the nature of transition to equilibrium is 
controlled by a parameter called by us ” Compton factor”, 
Kc = <?07, where 7 is a relativistic factor of electrons, 
0 = ksT/moc 2 is a normalized temperature, ks is the 
Boltzmann constant, rriQC 2 - the rest energy of electron, 
and q ~ 10, see below. When Kc C 1, we are deal¬ 
ing with a so called Thompson, or classical, scattering 
of light, whereby the scattering cross-section, <to is con¬ 
stant, even if 7 1, and the theory of light pressure is 

well known, see e. g. [5]. Due to advents in laser and con¬ 
trolled fusion technologies, the I<c can be large enough, 
Kc > 1, and we are entering a QED, or Compton do¬ 
main, where the respective theory is far less developed. 

In view of new developments in nuclear fusion and 
physics in general, e. g. in astrophysics and cosmol¬ 
ogy [29-31], it would be of fundamental importance to 
have the theory of that force over all the energies of 
electrons and temperatures of radiation, from the clas¬ 
sical, Kc -C 1, to transient, Kc ~ 1, to QED domains, 
Kc 1. From QM viewpoint [32,33], the light pres¬ 
sure on elementary particles is a result of averaging over 
ensemble of Inverse Compton Scattering events [34-36], 
whereby a particle transfers part of its momentum to a 
scattered photon. The known results provide patchy de¬ 
scriptions of the process, with quantitative results known 
mostly for “cold” case [5], and qualitative - for very “hot” 
case in the theory of high-energy cosmic rays [36-38]. Yet, 
to the best of our knowledge, no general formula for the 
light pressure F( 7,0) and related relaxation rates for ar¬ 
bitrary Kc presently exists. 

In this paper, 

(a) by using Lorentz transformation of spectrum, Sect. 
II, and relating photon scattering to momentum trans¬ 
fer to electrons via Doppler and Compton relationships, 
Sect. Ill, we derived a general light pressure formula, 
Sect. IV, for the isotropic/uniform radiation with an 
arbitary frequency spectrum, arbitaryly relativistic en¬ 
ergy of electrons and arbirary spectral dependence of the 
the scattering cross-section on frequency of light, 

(b) checked it out for the well known case of Thompson 
scattering and related light pressure, Sect. V, and sim¬ 
plified it in the case of blackbody radiation with Planks 
spectrum at arbitary temperature, Sect. VI; 
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(c) by using the frequency/energy dependence of 
photon-electron scattering cross-section in the Compton, 
or QED domain, based on the Klein-Nishina QED theory 
[39,40] which takes into consideration virtual electron- 
positron pair creation and annihilation, Sect. VII, we 
applied a general formula to the case of electron, im¬ 
mersed in the blackbody radiation, and found amazingly 
precise and universal analytic approximation for the light 
force in the entire domain, Sect. VIII, and finally 

(d) considered kinetics of electron density distribution 
to the equilibrium, i. e. the thermalization of the distri¬ 
bution, using Fokker-Planck equation and its solutions, 
in particular relaxation rates of the process, and how they 
may affect plasma oscillations, Sect. IX. 

Our results may have important applications for 
both high-temperature plasma, in particular controlled 
nuclear fusion and nuclear explosions, and to astro¬ 
physics/cosmology (to be addressed by us elsewhere [41]). 

II. LORENTZ TRANSFORMATION OF 
RADIATION SPECTRUM 

An isotropic (and homogeneous) radiation is associ¬ 
ated with some preferred L-frame, where a particle at 
rest experiences no time averaged light pressure, as the 
action of a fc-vector component is canceled by a counter- 
propagating (— k) component. We assume the spectrum 
of this radiation, p L (oj), known. A particle moves in the 
ai-axis in that frame with velocity v = ve x , and is at rest 
in a certain P-frame. A light pressure F = dp/dt = Fe x 
on a particle is then nonzero if v ^ 0 ; here p = pe x is its 
momentum in the L-frame, p/m^c = p = /Py = \f~f~- -T, 
where /? = v/c, m 0 is a particle rest mass, c - speed 
of light, and 7 = l/-\/l — /3 2 = \J\ + p 2 . Our deriva¬ 
tion of F(p) is based on Lorentz transformation of p{k) 
from L-frame to P-frame [42,43]; it is valid for arbitrary 
frequency dependence of a full cross-section <r(e) where 
e = huj/moc 2 , of scattering of an w-photon at a particle. 

For a gas of particles with non-zero mass, one can 
define distribution function g(p,r) in a lab L-frame in 
the phase space of momentum p and position vector r 
as the number of particles, dNpp = g(p,r)dfl , per the 
element of phase space, dfl = dVpdVp, where dVp = 
dp x dpydp z and dV? = dx dy dz are the elements of mo¬ 
mentum and coordinate spaces respectively. A general 
formula for a Lorentz transformation of a distribution 
function g L ( p L ,r L ) in the L-frame to a distribution func¬ 
tion g p ( p p , f p ) in a P-frame moving uniformly with re¬ 
spect to the L-frame is as [42,43]: 

= g P {p P ,r p y, dn L = dfi p (l) 

where k L and r L are related to k p and f p respectively by 
a standard Lorentz transform for an observer moving in 
L-frame in the a:-axis with velocity v = e x v. In the case 
of photon gas Eq. (1) remains true for the spectrum of 
photons g(k,r), by replacing p with k (k = to/ c), dVp by 


dVt = dk x dk y dk Zl and dfl by dfl = dV^dVp. Since we are 
interested here in the case of a homogeneous radiation, a 
spectrum is g(k), and its transformation is written as 

9l$ l ) = g P (k P ), dn L =dn p -, ( 2 ) 

(here g and dkl - dimensionless), where Lorentz trans¬ 
form for k is {k x ) L = 7 [{k x ) p + /3k p \, ( k±) L = ( k ±) p ; 
k L = 7 [k p + /3(k x ) p ], with fcj_ = k y e y +k z e z , k = w/c = 
\Jk x + k ' 2 ^, In particular, the Doppler coefficient is as 

D = uj p /uj l = [ 7(1 +/3cos £ P )] -1 = 7 (l-/ 3 cos£J (3) 

where £(...) are the angles between respective /c-vectors 
and the x-axis, i.e. cos (7 = ( k x ) L /k L , cos£ p = 
( k x ) P /k p , transformed as cos£ p = (cos (7 — /3)/(1 — 
/3cos£ p ). Furthermore, since the radiation is isotropic 
in the L-frame, the distribution function g L does not de¬ 
pend on the direction of k L -vector, and we have g L (k L ) = 
where k = \k\. Since both spectra are symmet¬ 
rical around the x-axis, we will use spherical coordinates 
in the k space, so that dVg = k 2 dk dO, where dO is the 
element of solid angle in the direction of k. One can then 
introduce the density number of photons in the element 
p(k), dk dOVp defined as 

p(k) = k 2 g(k) with dNj; = p(k)dk dO dVp; (4) 
([p] = cm -2 ), and its transformation using Eq. (2) as: 

P L ( k L)/ k l = Pp(kp)/k 2 p (5) 

From now on, since we are dealing with a homogeneous 
radiation independent on the radius-vector length, |r|, we 
will re-assign the notion of spectra p only to the ones be¬ 
ing functions of frequency u> and angles £ instead of k and 
r vectors. In this case, p(w,^)dw will have dimension of 
[cm -3 ]. As expected, the radiation becomes anisotropic 
in the P-frame (yet symmetric around the x-axis) with 
its spectrum p p (w, £ p ) given by a L-frame isotropic spec¬ 
trum p L (uj). whose argument and amplitude are now al¬ 
tered by the parameters j3 and £ p via Doppler coefficient 
D= [ 7(1 + (3 cos£ p )] -1 , Eq. (3): 

P P {w^ p ) =p L [u/D((3,t P )]D 2 (p,t; p ) ( 6 ) 

Our further calculations will mostly be focussed on the 
events in the P-frame assuming that p L (w) is known. 

III. PHOTON SCATTERING AND 
MOMENTUM TRANSFER 

We designate the wave-vector of an incident photon 
in the P-frame as ki n , and that of a scattered photon 
as k sc . The latter one is scattered into the solid angle 
dO sc = sin il) sc dil> sc d(j) around k sc . Here ij) sc £ [0, 7 r) is 
the angle between incident and scattered fc-vectors, and 
(j) £ [ 0 , 27t) is an azimuthal angle around the k sc direction. 
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The number of photons scattered into dO sc within time 
interval dt and spectral band dto = cdk , is 

dN sc = p p (uJi n , £in){dcr/ dO S c)(dO SC /4n)cdtdto (7) 


where a MT is the cross-section of a momentum transfer 
from Wj„-photons to a particle. It must be noted that 
a MT may not in general coincide with a plain f ull (inte¬ 
grated) scattering cross-section, a(ci n ) 


where the differential scattering cross-section da/dO sc 
describes an (unknown yet) physics of energy and mo¬ 
mentum transfer from the incident photon to both scat¬ 
tered photon and a particle [including a possible excita¬ 
tion of internal degrees of freedom in the particle if there 
is any, which is not the case for an electron, whereby 
da/dOgc is strictly due to Compton scattering, see next 
equation ( 8 ), that enters Klein-Nishina formula, see be¬ 
low, Eq. (27), and discussion in the end of Sect. IV]. 

The Compton quantum formula determines the ratio 
R of photon energies after and before scattering from a 
single electron: 

Csc/ejn = R{ein,ipsc) = [1 + e»„(l - cos Vtsc )] -1 ( 8 ) 


a(e in ) = 27r 



da {tin) 
dO sc 


sin ip sc dips C ; 


tin 


mgC 2 


(13) 


because only the projection of k sc into fcj ra -direction con¬ 
tributes to the radiation force F, see Eqs. (9) and (11), 
where in general R(e,ip sc ) < 1 if e > 0, Eq. ( 8 ). a MT 
and cr(e,; n ) are related as a MT = a — an, where 


/»7r do" 

a R (e) = 2ir / R{e,ip)——cosip sin ip dip = 
Jo dO 


27r 



da pdp 

dO 1 + e(l - C) 


(14) 


where e = hoo/moc 2 . When calculated back to the L- 
frame, photons back-scattered from an ultra-relativistic 
electron, 7 1 , may have large energies with the 

Doppler shift up to D « 2y, even if their energies in 
the P-frame are still below QED limit, i. e. Cm <C 1. It 
is commonly called an Inverse Compton (or Thompson, 
if -C 1 ) scattering. 

Only the projection of k sc into /^-direction contribute 
to the force F ; all the rest are canceled out after integra¬ 
tion over the azimuthal angle </> (in the P-frame, where 
the electron is at rest, the scattering problem has a sym¬ 
metry around kin). Thus the momentum transfer to an 
electron in /^-direction after the scattering is 

hAktr = Hhn-ksc cos ip sc ) = htOi n [l-R(e,ip sc ) cosip sc ]/c 

(9) 

IV. LIGHT PRESSURE ON A PARTICLE 


Considering the number of photons, dN sc , Eq. (7), 
scattered into a solid angle dO sc , the light pressure im¬ 
pacted by them on the electron in fc,; n -direction, is the 
rate of momentum transfer, F = Fcfin where qi n = 
hn/hn and 


dF = hAktr 


d 2 N s , 

dtdu>i 


hujin da 

_^ 

47T dOsc 


[1 - R(e, ipsc) cos ip sc }p p (u)in, fin) sinipscdipscdfp (10) 

Integrating Eq. (10) over <p and ip sc , we find a full u>. in- 
Fourier component of the light pressure as: 

H ^in-> Pi Qin ) — [hiJin/ 47t) p p (u7 n , fin)a MT (H) 


with 


r 

a MT =27r / [1 - P(e, ip sc ) cos ip sc ] 777 - sin ip sc dip s 
Jo dO sc 


( 12 ) 


which in turn reflects P-factor in Eq. (9) and (12). It 
zeroes out for a = const (Thompson scattering) and at 
e = 0, and peaks at e ~ 0.54 (see the end of Sect. VII), 
but is negligibly small both at e <C 1 and e 1. 

Now, we compute the light pressure, P = F{p)e x , as 
an integral of the ^-projections of the Fourier force com¬ 
ponents, F(u)in,P,q in ), i- e. T x = F(to in , /?, q in )e x = 
F cos fin , over all the incident solid angles, dOi n = 
sin fi n dfindcpin, and frequencies u>i n in the P-frame. 
Thus, we have for the full light pressure 



™ ■ t - 


dF x 

dOir, 


■dOindlO = 


F , pi , Cm ) SIC fin COS findfii 


dto 


(15) 


Recalling that a spectrum p p in F, Eqs. (11),(15) can 
be expressed directly via the known isotropic spectrum in 
the L- frame, p L , Eq. ( 6 ), we can now write the expression 
for the light pressure in closed form as: 


h f°° 

F(p)=~l WCr MT (w) 

(16) 

or by using a substitution, f = cosf , reduce it to: 

„ n = ^ r _ <Jmt ^ I f 1 Pl\fi { 1 + PQ\ 

W 2 Jo 7 2 U-i (1 + /3C) 2 

(17) 

Alternatively, by using a substitute v = wy(l + pip), the 
same result can be written as 


pdf j du) 


sin(2 f)p L ( u)/D)D 2 df 


dto 


(18) 

Note also that only full (integrated) cross-section, a MT , 
enters into the final calculations. Which one of Eqs. (17) 
or (18) to use for detailed study is a matter of computa¬ 
tional convenience depending on specific model functions 
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v h 
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it mt (s) and p L (v) (see Sect. VI below for Planck radia¬ 
tion) ; both of them incorporate ensemble averaging over 
all the relevant parameters, which makes pressure F the 
best tool to explore electron de-acceleration in EM field. 

Let us reflect on the domain of validity of Eqs. (17) 
and (18). The rest of this paper is dealing with the 
light pressure on a single, most fundamental elementary 
particle, electron (or substantially rarefied electron gas 
or plasma), which presents a clear case. The question 
is then whether they could be applicable for more gen¬ 
eral cases, in particular for high density gas or plasma 
with many-body interactions, in particular plasma oscil¬ 
lations, and for single particles/objects with an internal 
structure and resonances. We address the former issue in 
the Sect. IX below [see the text preceding and following 
Eq. (39)], and discuss the latter one here. 

In our derivation of Eqs. (17), (18) we have not used 
any assumption based on the fact that a scattering ob¬ 
ject is an elementary particle. Essentially, the only as¬ 
sumption was that we have only one particle and one 
photon both in input and final output channels in each 
act of scattering, and were not concerned about interme¬ 
diate processes. The physics of these processes in gen¬ 
eral is to be described by the differential cross-section 
d<j(ui)/dOgci which in the case of an electron is due to 
Klein-Nislrina formula, see below Eq. (27), based on 
Compton scattering relationship for input/output elec¬ 
tron energies, Eq. (8); eventually, da/dO sc is absorbed 
into the full cross-sections of scattering er(w) and momen¬ 
tum transfer, a MT (w) via integration, Eq. (13) and Eq. 
(12), which may also include a new ratio R of photon 
energies after and before scattering from a particle, that 
may differ now from the Compton quantum formula, Eq. 
(8), for a single electron, based now on radiation loss or 
amplifications after scattering from the particle due to 
its internal degrees of freedom. Thus, Eqs. (17), (18) 
use only our knowledge of a and p as functions of ui. 
They will be valid for example for the resonant or other 
dispersion-related interaction of the radiation with atoms 
or even macro-particles, with any quantum or classical 
resonances due to e. g. dipole momenta, band struc¬ 
ture, eigen-modes, etc. By the same token it also does 
not matter whether the cross-section a is due to elastic 
scattering or includes losses of energy to internal degrees 
of freedom; all that implicitly enters into the functions 
cr(cv) and cr MT (ui). To a degree, this is reminiscent of 
a phenomenological role played by dispersive dielectric 
constant e{ui) in electrodynamics whereby e(u>) provides 
a short-hand representation of all the constitutive inter¬ 
actions of EM-field with matter. 

The above discussion clearly suggests the situations 
whereby the description offered by Eqs. (17), (18) is not 
comprehensive: it is when there are more than one output 
particle (as e. g. in the case of photoionization resulting 
in an ion and one or more ionized electrons), and/or more 
than one output photon (as e. g. in stimulated emission 
or any kind of nonlinear multi-photon process, such as 
sum, difference, or high harmonics generation, etc). In 


all those cases, one needs to include all the scattering 
channels and generalize Eqs. (17), (18) by summation 
over all of them. 

V. THOMPSON LIGHT PRESSURE 

(°mt =cr = const, a R = 0 ) 

In the limit of frequency-independent cross-section, the 
integral in Eq. (17) is readily evaluated resulting in: 

F Th = -(4/3)ct 0 W l /z7, with 7 = \/l + p 2 (19) 

where W L = h / 0 °° u>p L (oj)du) is the energy density of 
photons in the L-frame, and ao is a full (classical) cross- 
section of a charged particle: 

a = const = <7 0 = (87t/ 3)r 2 ; with r 0 = e 2 /m 0 c 2 (20) 

where r 0 is a classical EM-radius of a particle. Eq. (19) 
coincides with known results [5]. At p <C 1 we have 
—F oc p, while in a relativistic case, p 1, — F oc p\p\. 
It is reminiscent of a drag force in liquids and gases, which 
is linear in velocity v for low v (Stokes force), and oc v|v| 
for highly turbulent flow [44]. 

VI. LIGHT PRESSURE OF A BLACKBODY 
(PLANCK) RADIATION 

Spectral, p TR , and total energy, W TR = h / 0 °° wp TR (w) 
du> densities of blackbody radiation in the L -frame at the 
temperature T (« 2.725K for current CMB) are 

, x , W 2 dw 8n 5 4 

p TR {w)du — W TR — ^ W c 9 

.' , (21) 
which is a familiar Planck density distribution, where k B 

is the Boltzmann constant, W c = m 0 c 2 /X^, is a “Comp¬ 
ton energy density”, A c = 2Trh/moc is the Compton 
wavelength, and 9 = k B T/m 0 c 2 is a dimensionless tem¬ 
perature (for the current CMB, 9 « 0.534x 10 -9 ). In 
Thompson limit, the energy density W L can now be re¬ 
placed by W TR . For further calculations, we will use a 
dimensionless time r = t/t c , and force / = Ft c /n%oc by 
introducing a “Compton time scale” for an electron: 

t c =135A c /647r 4 a 2 c«3.25x 1CT 18 s; t c oc h 3 (22) 

where a = e 2 /he « 1/137 is the fine structure constant. 
(It is worth noting that a “t/-scale” 9 V = {t c /t u ) 1 / i ~ 
1.65 x 10~ 9 , where t v ~ 4.4 x 10 17 s is the age of the 
universe, comes close to the current CMB temperature, 
9 « 0.534 x 10~ 9 .) In dimensionless terms, Eq. (19) for 
Thompson limit can now be rewritten as: 

(■ dp/dr) Th = f Th = -9 4 p<y (23) 

This approximation is valid for Kq 1 (hence 9 <C 1), 
and it is still good for relativistic case, p ~ 7 » 1 , as long 
as p <C 9~ l . (Note that for electrons, 9 = 1 corresponds 
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to T ps 0.6 x 10 lo .A). Eq (23) is readily solved for /x(r); 
with an initial condition fi = /x 0 at r = 0, we have 


1 , 4 4 

2 + 7 ~ 2(1+ 2e) 2 


. 7 hu 

with e =-- 

m Q c z 


(28) 


H(t ) = l/sinh(rd 4 + <5 0 ); (24) 


where <5 0 = Zn[(l+7o)//xo]. At /x <C 1 Eq. (24) reduces to 
/x oc exp (—t/t TR ), while in relativistic case, fi 1, - to 
H « /j, 0 (t/t TR + 1) _1 ; a time scale here is t TR = t c / 7 0 d 4 . 
This scale may vary tremendously even for /xo < 1 - 
from 0.4 x 10 _10 s = 40ps for the temperature ~ 10 8 K 
(9 ~ 1.7 x 10“ 2 ) below lab-nuclear fusion, to 0.4 x 10 2O s 
- for the current epoch CMB (9 ss 0.534 x 10 -9 ), which 
is 10 2 times longer than the age of the universe, t v . 

For a frequency-dependent cross-section <j MT (e) , where 
e = Uuj/moc 2 , Eq. (17) for a blackbody radiation (21), 
can be readily reduced to a single integral by using the 
dilogarithm function, L^(z) = — / 0 “ ln( 1 — t)dt/t [45]: 


/ 


45 6» 4 

87 t 4 (/xy) 2 



°MT M/t) 


(5“ + S + )dx-, 


(25) 


where 

S± 




e -s(l±/3) 


x 2 fd In 


2 _ e -x(i±/3) 


and f3 = /x/y. In particular, for a low-relativistic motion, 
(3 ~ /x <C 1, but arbitrary high temperature 0, Eq. (25) 
is further reduced to f = — /x0 4 0(0), where 


0(d) 


/"” <j mt {x6)x 4 
167t 4 J q tr 0 sinh 2 {x/ 2) 


0 ( 0 ) 


1 (26) 


For arbitrary yx, a specific case of electron is considered in 
Sect. VIII below, but we need first to determine a QED- 
related energy dependence of scattering & momentum- 
transfer cross-sections for electron in the next Section. 


VII. QED SCATTERING & 
MOMENTUM-TRANSFER CROSS-SECTIONS 
FOR ELECTRON 


In the “cold” and “hot” limits we have respectively, 
<x KN /a 0 « l-2eat e < 1; and a KN /cr a « (3/8e)[ln(2e)+ 
1/2] at e^> 1. Using now Eq. (8) in Eq. (14), we have 

— = lf 1 &[& - (1 - C 2 W + l]CdC (29) 

CT 0 a 7-1 


where R = [1 + e(l — £)] 1 - Its integration yields: 


CT„(e) 3 f 2e(l — 2e) 

2(1 + e) ln(l + 2e) 

a 0 8(3(l + 2e) 3 

e(l + 2e) e 2 


f 2(2 + 4e + e 2 ) 3(1 -I- e) 2(1 + e) ln(l -|- 2e) 

\ e(l + 2e) 2 e 3 [ 1 + 2e e 

(30) 

As mentioned already, cr R (e = 0) = 0. The three terms 
within the “outer” brackets {...} in Eq. (30) are grouped 
to have each one of them to also zero out at e = 0. In the 
“cold” and “hot” limits we have respectively a R /a 0 « 
6e/5 at e <C 1, and cr R /cr 0 « (l/2e)[l — 31n(l + 2e)/4e] at 
e 1, so that the term a R in the photonf->electron mo¬ 
mentum transfer can be neglected at e,£ _1 <C 1. It can 
be shown that a KN > a R > 0 everywhere in e £ (0,oo), 
so that total cross-section a MT = cr KN — a R in Eqs. (17) 
and (18) is always positively defined. All the ct’s spectral 
profiles are depicted at Fig. 1 showing that a R (e) peaks 
as a R /a 0 « 0.1442 at e w 0.543 {Twj « 277.5 KeV). 
At that point a KN « 0.55u o , and a MT ss 0.4<r o , i. e. 
ct k /&kn ~ 0-26 a. Thus, strictly speaking, cr R should 
not be neglected within KN-theory using O(a) terms, at 
least around e ~ 1. Yet in reality, it makes little dif¬ 
ference when calculating /(/x, 9) in the entire momentum 
span, /x € (0,/x Pi ), where [i pl = ksTpi/rrioc 2 « 2.4 x 10 22 
is the highest momentum in the universe related to the 
Planck temperature Tpi ~ 1.417 x 10 32 A', see below. 



In the limit of a low-energy photons, e <1, their scat¬ 
tering by a charged particle is described by an energy- 
independent Thompson cross-section ao, Eq. (20). 
Yet a cross-section a becomes energy-dependent even 
at sub-relativistic energies, which in the case of elec¬ 
trons/leptons is due to quantum Compton scattering 
via Klein-Nishina theory providing an exact solution for 
a = a KN for any e, good to the first degree in a. The 
differential cross-section da KN /dO in that case is [39,40]: 

da KN /dO = (3a 0 /16n)R 2 (R + A -1 - sin 2 (27) 

with R = R{e,ip) as in Eq. (8). Using Eq. (27) in Eq. 
(13) we get a full Klein-Nishina cross-section [39,40]: 


CT™(e) 


G 


o 




ln(l + 2e)+ 


VIII. QED BLACKBODY LIGHT PRESSURE 

Eqs. (28),(30) together with Eq. (24)-(26) allow for 
specific investigation of the light pressure by Planck ra¬ 
diation on an electron. The force / vs /x in the entire 
momentum span, /x £ (0,/x pi ), and for various 8, from 
9 ~ 10 -9 to 10 6 , was numerically evaluated using Eq. 
(25) and depicted in Fig. 2 for the relativistic/QED fac¬ 
tor Q(n,Q) = — f/n9 A , where — /xd 4 is the Thompson 
non-relativistic light pressure, Eq. (23), for 7 = 1; for 
/x, 0 <C 1, Q = 1. In ultra-relativistic QED, or Compton 
domain, based on the behavior of a KN (w) at ui —>• 00, Eq. 
(28), the asymptotics of Q at /x > maa;(l,0 -1 ) can be 
shown to be Q oc [ln(0/x) + 0(\)\/9. Using these numeri¬ 
cal and asymptotic results, it would be greatly beneficial 
for further analysis to have their good analytic interpo¬ 
lation/approximation. Amazingly, this task is perfectly 
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FIG. 1: Dimensionless integrated cross-sections of photons 
scattering by an electron: Klein-Nishina a KN /a 0 , “pro¬ 
jection” term cr R /a 0 , and momentum-transfer a MT /a 0 = 
( a KN ~ <j r)/ <j oi vs dimensionless energy of incident photons. 



FIG. 2: QED/relativistic factor, Q = \f\/p9 4 vs momentum, 
p, for various temperatures, 9. Three curves slightly divergent 
around p ~ p c in the case 9 = 9 CMB correspond to numeric 
integration of Eq. (17) with a MT = a KN —cr R (upper), a MT « 
a KN (middle), and analytic Eq. (32) (lower). Similar curves 
for 9 = 10 6 coincide to the line width. Inset: a QED factor 
0, Eq. (27), for non-relativistic motion vs 9. 

served by a remarkably simple formula good for the entire 
span p e (0, hpi) and 9 e (0, 0 Pl ): 


IX. KINETICS OF DENSITY DISTRIBUTION 
NEAR TO AND FAR FROM EQUILIBRIUM 

Since F = dp/dt , Eq. (15), or / = dp/dr [see e. g. Eq. 
(23)], the averaged dynamics of electron motion, p(r), for 
9 = const is implicitly described as r = f dp/ f(p,9). It 
is readily integrated in the case of Thompson scatter¬ 
ing, Kc <C 1 yielding an explicit function p(r), Eq. (24), 
whereas in general case, especially for the transition from 
Compton domain, Kc 1 to the Thompson domain, the 
timeline becomes more complicated, yet still analytically 
solvable using another approximate model, which is very 
close to the Eq. (31) in the energy span covering the en¬ 
tire Compton and most of the Thompson domains [41]. 
In the context of this paper, an important issue is the re¬ 
laxation time T r i x (or rate t~ 1x ) vs temperature T of the 
electron distribution to its equilibrium state at a given T, 
and the momentum p of a non-equilibrium electron. This 
problem is best handled by using a Fokker-Planck equa¬ 
tion for the diffusion in the momentum space [46]. To 
that end, we consider a distribution function, g^ e \p,r) 
of electrons defined here as the number of electrons per 
element of solid angle dO , element of momentum, dp, 
within a unity of coordinate space, and a density num¬ 
ber, // e )(/x,T) = 4irp 2 g( e \ Note that in application to 
cosmology, whereby one needs to consider the expanding 
space/universe, these functions reflect the distribution 
within the expanding unity of coordinate space, i. e. in 
the spatial “unity box” expanding at the same rate as 
the Universe. Assuming then that (a) the electron distri¬ 
bution is isotropic, same as CMB, (b) the total number 
of electrons, in the unity of momentum space and the 
expanding unity of coordinate space is approximately in¬ 
variant, / 0 °° dp = iVs = inv, and (c) the thermal 
equilibrium of a relativistic gas at any 9 is due to the 
Maxwell-Jiittner (MJ) distribution [47,48], 

g^ j <xe-^ e [9K 2 (l/9)}- 1 (33) 

where K 2 is the modified Bessel function of the second 
order, with MJ being a relativistic generalization of the 
Maxwell-Boltzmann distribution, 

9^ B oce-^ l26 9- 3 / 2 (34) 


Im(p, 9) = -p9 3 ln (1 + K c )/q\ K c = 7 9q (31) 

where q = 10 is a fitting parameter. Eq. (31) makes 
better than a few percents fit to the numerics over the 
entire span of momentum but a small area near Kc ~ 1 , 
Fig. 2, where they are still very close; it can be viewed 
as a benchmark for any other possible approximations. 
In the limit Kc 1, Eq. (31) is reduced to Eq. (23). 
For low-relativistic motion, p <C 1, 7 « 1, yet arbitrary 
temperature, the factor 0(0) in Eq. (26), shown at the 
inset in Fig. 2, is now approximated using Eq. (31) as 

0 = ln( 1 + 9q)/9q; (0 ss 1 at 9q <C 1) (32) 


we found a Fokker-Planck equation for the distribution 
function g( e \p,r), in terms of dimensionless momentum 
p, factor 7 , time r, temperature 9, and force /, as 


^V e) ] , d , 2 „ , 

- g — + -PfC,r) 


9 V 


7 <9o (e) 
0(r) 2 ^- 
11 0(1 


= 0 


(35) 

This equation is valid for time-dependant 9 and for the 
expanding universe For our purposes here we will con¬ 
sider only the case 9 = const = 9 eq (and neglect the 
universe expansion), since the time involved is mush 
shorter than the age of universe for the effects of inter¬ 
est here. Eq. (35) is solved fully analytically for the 
non-relativistic case, p <C 1, whereby any initial distri¬ 
bution function g^ (r = 0) is decomposed into Gaussian 
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components, Eq. (34), each of which has its specific ini¬ 
tial temperature, 9^ at r = 0; in time, their distribution 
functions remain Gaussian, with their time-dependant ef¬ 
fective temperature being as: 

0 (e) (r) = 9 eq + (*£> - 0 eq )e ~^, TrJx = 1/2 9% (36) 

i. e. t r i x = t T R/2 = t c /29j q (see eq. (24) and ex¬ 
planations therein). Eq. (36) still holds approximately 
for the entire Thompson domain for initial conditions 
close to the equilibrium (be reminded that the condition 
Kq 1 automatically means 9 eq <C 1, while proximity 
to the equilibrium - that 7= 0(1) ). In general case of 
arbitrary initial conditions the relaxation rate becomes 

T rll ~ 9 _1 [26'eq^( 1 + ^eqlin)\\ (37) 

where 7is the energy at the peak of initial density 

(e) 

distribution. If it was MJ density distribution p M j oc 
°f the initial temperature 9i n , then 74„ = [1 + 
29 in {9 in - l-i/l + 9j n )] 1/2 , and t rix = t C T r i x . For 9 in < 1, 
it coincides with Eq. (36), and for 9 2 n > 1 we have 

Tri x » q[29l q ln{2q9 eq 9 m )]- 1 . (38) 

It is instructive to look at a few examples of interest. 
For the current CMB, the relaxation time to the equi¬ 
librium will exceed the age of universe by ~ 10 2 , see. 
Eq. (23) and related discussion, which makes the ther- 
malization here almost irrelevant. With an easily lab- 
accessible T = 1300A' (9 = 2.2 x 10~ 7 ), we have t r i x » 
0.7 x 10 9 s, which is still unrealistic in practical terms. 
At the ’’Compton threshold”, Kc = 1, 9 = 0.1, we have 
t rlx « 2.3 x 10 _14 s = 23/s. With T ~ 10 lo A' (9 « 1.7) 
required for the controlled nuclear fusion [49,50], we have 
trix ~ 10~ 18 s = las for 7~ 1. If 7j n ^ 1, the pro¬ 
cess is getting even faster, and the radiation might act to 
almost instantly deplete coherency of e. g. an electron 
beam used as a plasma diagnostic tool. For example, if 
its energy is 50MeV, a Compton factor in Eq. (37) is 
Kc = q9 e q r )in ~ 1,700, with t r i x ~ 0.4as, so that the 
starting rate of e-beam thermalization is much faster. 

Having in mind plasma nuclear fusion, it would be of 
substantial interest to see how the light pressure-induced 
damping may affect plasma oscillations in high-density, 
high-temperature plasmas. Those many-body excita¬ 
tions seem to preclude a single-electron light pressure, 
Eq. (31), from playing a significant role in plasma re¬ 
laxation. However, the relaxation time t r i x due to light 
pressure does become a major player as long as it gets 
faster than relaxation in other channels of thermalization 
such as collision with similar or other species (e. g. elec¬ 
trons and protons), t c i s , to dominate in a total relaxation 
rate t-^ t r i x T t c ^ s if trix t c i s , hence t^ ~ t r i x • Yet 
the most important effect here is that the time t r i x may 
get even much shorter than a plasma oscillation period, 
trixLOpi < 1, in which case those oscillations would be 
damped or even extinguished. Using in rough approxi¬ 
mation a standard equation for plasma frequency oj p i vs 


the number density of electrons, N e , we get the condition 
on N e to have plasma oscillations strongly damped 


N e < N cr 


mo/ire 2 
( t'c Trix) 2, 


1.2 x 10 26 


-cm 


-3 


‘ rlx 


(39) 


up to w ~ Udmp = trix’ where K cr is a critical number 
density of electrons below which the light pressure sup¬ 
presses plasma oscillations. For T ~ 10 lo A' (9 = 1.7) 
required for the nuclear fusion [49,50], presuming near¬ 
equilibrium, 9 eq S3 9i n , and using Compton domain Eq. 
(38), we have N cr ss 2 x 10 27 cm -3 , which exceeds the 
number density of the interior of most of the stars [29- 
31], the Sun incuding, thus making plasma oscillations 
of electrons completely extinguished. Even with an or¬ 
der of magnitude lower temperature T ~ 10 9 A', using 
Thompson formula (36) for r r i x in Eq. (39), we get 
N cr S3 3.4 x 10 2 °cto - 3 , which is still much higher than 
any conceivable lab plasma density. 

Notice that far from equilibrium, with 9i n 9 eq in 
Eqs. (37), (38), the thermalization rate, as well as 
the ’’friction coefficient”, f/p in Eq. (31), still increase 
with energy 7j„ (which translates into photon energy e in¬ 
crease in P-frame) in the Compton domain, Kc 1, for 
a fixed temperature, while cross-section a(e) decreases. 
The explanation of this is that while the <r(e) is indeed 
slowly receding with photon energy (as ~ l/ln(e), Eqs. 
(28) and (30)), each act of Inverse Compton Scattering 
(ICS) gets much more “quantum efficient” since then a 
low-energy photon scattered from a high-energy electron 
gets a huge boost by accruing up to almost full energy 
of the electron. The peak gain is reached in a “head-on” 
collision, when a photon is exactly back-scattered. Based 
on the Compton scattering formula, Eq. (8), and Doppler 
effect, Eq. (3), that enter the final formula, Eq. (17), the 
maximum scattered photon energy in the L-frame is 

e sc = Tiu sc /moc 2 = e in (7 + p)/(j - p + 2e in ) (40) 


where <7 n is the incident photon energy. For high-energy 
electrons, 7 ss p 1, and low-energy incident photons, 
tin <C 1, the maximum quantum efficiency of ICS de¬ 
fined as the ratio of scattered photon energy to that 
of an incident electron, r) = e sc /7 , is then rj « [1+ 
(47 d%n) ] 1 and in the sub QED domain, 7 1/ din, 

we have 77 « 47c<C 1, hence small loss of electron en¬ 
ergy per collision. However, in QED Compton domain, 
4y din 1, we have 77 « 1 — (47ei n ) -1 ~ 1, i. e. an elec¬ 
tron passes great part of its energy to a scattered photon. 
With considerable probability an electron jumps then in 
one collision to the Compton threshold, p —t p sc ~ 1 /9q. 

These effects make it of special interest to look into the 
kinetics of highly relativistic electrons with their energy 
far exceeding that of equilibrium. Since in such a case 
the system remains far from the equilibrium during the 
evolution, the last term in Eq. (35) can be omitted, and 
in terms of density p^ oc p 2 g^ it can be reduced to 

dp {e) /dT + d[fp (e) }/dp = 0 


(41) 



which is essentially a continuity-like equation and is fully 
integrable; its general solution can be shown to be 

P (e) = $(£ - with £ = f dfi/f (42) 

where < f>(s) is an arbitrary function of s ; here it is de¬ 
fined by initial conditions, e. g. the MJ-distribution 
with t > 1. Resulting evolution of the spectra of high- 
T sources in the relict radiation environment reveals their 
drastic transformation, e. g. formation of narrow spec¬ 
tral lines in the cosmic electron spectrum, and develop¬ 
ment of a ” frozen non-equilibrium” state due to very low 
rate of electron momentum decay in the Thompson do¬ 
main. The ramifications of these effects in astrophysics 
and cosmology will be discussed by us elsewhere [41]. 

X. CONCLUSION 

In conclusion, we developed a theory for the light pres¬ 
sure on particles, in particular electrons, by isotropic ra¬ 


diation, in particular blackbody/Planck radiation, that 
covers the entire span of energies/momenta up to the 
Planck energy by using, in the case of electron, the QED 
Klein-Nishina theory for electron-photon cross-section 
scattering. We also analyzed the kinetics of electron 
relaxation into equilibrium by using relativistic Fokker- 
Planck equation for temporal evolution of electron spec¬ 
tra of high-T sources, revealing dramatic difference be¬ 
tween classical and QED electron relaxation rates that 
may result in a host of new effects. We showed as well 
that a light pressure-induced damping may completely 
extinguish plasma oscillations of electrons at the tem¬ 
peratures approaching nuclear fusion. 
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